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The pair-breaking critical current density of magnesium diboride 
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The pair-breaking critical current density, jd, of magnesium diboride was determined over its 
entire temperature range by a pulsed dc transport measurement. At fixed low values of current 
density j, the resistive transition temperature T c shifts in the classic ATc(j)/Tc(0) tx — \j/jd(0)] 2 ^ 3 
manner, with a projected jd{0) m 2 x 10 7 A/cm 2 . The directly measured jd(0), from current- 
voltage (I-V) curves at different fixed temperatures, has a similar value and the overall temperature 
dependence jd(T) and magnitude are consistent with Ginzburg-Landau theory. 

PACS numbers: 74.25.Sv, 74.25.Fy, 74.25.Bt 



INTRODUCTION 

Magnesium diboride (MgB2) recently made an impact 
as a promising new superconductor with a surprisingly 
high critical temperature for a simple binary compound. 
This has spurred considerable research activity into in- 
vestigating the myriad properties associated with its su- 
perconducting state. Besides the critical temperature T c 
and the upper critical field H C 2, an intrinsic parameter 
that sets a fundamental limit to the survival of supercon- 
ductivity is the pair-breaking (or depairing) critical cur- 
rent density jd- We report the first measurement of this 
important quantity in the MgB2 superconductor, which 
sets an absolute limit to the maximum current-carrying 
performance under ideal conditions. This also represents, 
to our knowledge, the only complete (0 < T < T c ) mea- 
surement of jd by a direct transport method in any type- 
II superconductor. 

When a superconducting state is formed, charge car- 
riers correlate and condense into a coherent macroscopic 
quantum state. The formation of this state is governed 
principally by a competition between four energies: con- 
densation, magnetic-field expulsion, thermal, and kinetic. 
The order parameter A, that describes the extent of con- 
densation and the strength of the superconducting state, 
is reduced as the temperature T, magntetic field H, and 
electric current density j are increased. The bound- 
ary in the T-H-j phase space that separates the super- 
conducting and normal states is where A vanishes, and 
the three parameters attain their critical values T c (H,j), 
H c2 (T,j), and j d (T, H). 

In practice, a superconductor loses its ability to carry 
resistanceless current long before j reaches jd- Any pro- 
cess that causes the phase difference between two points 
to change with time — such as the motion of flux vor- 
tices, phase slip centers in narrow wires, junctions, and 
fluctuations — can generate a finite voltage and hence re- 
sistance. The conventional critical-current density j c 



marks this onset of dissipation — depending on extrin- 
sic variables such as vortex pinning by defects — and in 
type-II superconductors can be a few orders of magni- 
tude lower than jd- Thus the transport becomes resistive 
and intensely dissipative long before the thermodynamic 
limit is reached, tending to mask a direct measurement 
of jd by sample heating. In this work we use a highly 
evolved pulsed-current technique (which we have refined 
over ten years) to reduce heating and obtain jd{T — > 0) 
both from (1) a direct measurement of j required to drive 
the system normal at T <CT C and (2) from the shift in T c 
as a function of j near T ^T c . 



THEORETICAL BACKGROUND 

A theoretical estimate of jd can be obtained from the 
Ginzburg-Landau (GL) theory, in which the strength of 
the superconducting state is expressed through the com- 
plex phenomenogical order parameter -0 = \ip\e lip . The 
superfluid density near T c is proportional to \ip\ 2 and 
the free-energy density / of the system (w.r.t. the free- 
energy density in the normal state) can be expressed as 
a power expansion in \ijj\ 2 (In "dirty" superconductors — 
superconductors with a high impurity scattering rate — 
the approximate validity of the GL expressions extends 
down to T <C T c .). In the absence of significant magnetic 
fields and in situations where the magnitude of the order 
parameter \ip\ is uniform (either because the dimensions 
of the sample are small compared to the coherence length 
or because of the principle of minimum entropy produc- 
tion at high dissipation levels Q) / can be expressed asQ 



/ = a|^| a + ||V| 4 + i|V| a mX a - 



(1) 



a and (3 are negative and positive constants respectively 
(a becomes positive above T c ), and the positive third 
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term is the kinetic energy density expressed in terms of 
the superfluid velocity v s — ^J£. — where e* and 

to* are respectively the effective charge and mass of a 
cooper pair. For zero v s , the equilibrium value of \xjj\ 2 
that minimizes the free energy (Eq. |l|) is IV'ool 2 = —a/ (3. 
For a finite v s it becomes 
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The corresponding supercurrent density is 
j = e*\i>\ 2 v s =2e|?Aco| 2 (1 
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The maximum possible value of this expression can now 
be identified with j d : 



j d {T) = 2e|^c 



2a 

3 V3^* 



1/2 



cH c (T) 
3V6ttA(T) 



(4) 



where the GL-theory parameters were replaced by 
their expressions a(T) = -(e* 2 /m*c 2 )H 2 (T)X 2 {T) and 
(3(T) = (47re* 4 /m* 2 c 4 )ff c 2 (T)A 4 (T) in terms of the phys- 
ically measureable quantities H c (thermodynamic crit- 
ical field) and A (magnetic penetration depth). The 
relatio ns H C (T) » H C (0)[1 - (T/T c ) 2 } and A(T) « 

A(0)/vTl - WJW\ g^e 

« irf(0)[l ~ (T/T c ) 2 ]t[l + (T/T c ) 2 }^ (5) 

where 

j d (0) =ci7 c (0)/[3V6vrA(0)] (6) 

is the zero-temperature depairing current density. Close 
to T c , Eq. | reduces to Jd (T » T c ) a 4j d (0)[l - T/T c ]i. 
This can be inverted to give the shift in transition tem- 
perature T C (J) at small currents, with the well-known 
j2/3 p r0 p OT tionality: 

T c (0) - T c (j) 
Tc(0) 





j 




Ud(o)\ 



(7) 



(The preceding discussion is based on Refs. || and ||.) 
Note that if heat removal from the sample is ineffective, 
Joule heating will give an apparent shift AT C oc pj 2 , 
which is the cube of the intrinsic ~ j 2 / 3 depairing shift 
near T c , and hence easily distinguishable. 



EXPERIMENTAL DETAILS 

The samples are 400 nm thick films of MgB2 fabri- 
cated using a two-step method whose details are de- 
scribed elsewhere^, |^]. An amorphous boron film was 
deposited on a (1102) AI2O3 substrate at room temper- 
ature by pulsed-laser ablation. The boron film was then 
put into a Nb tube with high-purity Mg metal (99.9%) 



and the Nb tube was then sealed using an arc furnace 
in an argon atmosphere. Finally, the heat treatment was 
carried out at 900° C for 30 min. in an evacuated quartz 
ampoule sealed under high vacuum. X-ray diffraction 
indicates a highly c-axis-oriented crystal structure nor- 
mal to the substrate surface with no impurity phases. 
The films were photolithographically patterned down to 
narrow bridges. In this paper we show data on three 
bridges, labelled S, M, and L (for small, medium, and 
large) with lateral dimensions 2.8 x 33, 3.0 x 61, and 
9.7 x 172 /1m 2 respectively. The lateral dimensions are 
uncertain by ±0.7/xm and the thickness by ±50 nm. 

The electrical transport measurements were made us- 
ing a pulsed signal source with pulse durations ranging 
0.1-4 /zs and a duty cycle of about 1 ppm. From past ex- 
perience with other films (e.g., YiB^CuaOy on LaAlOa) 
we found that micron-wide bridges typically have thermal 
resistances of order R t h ~ 1-10 nK.cm 3 /W at microsec- 
ond timescales [fil H, mj. For the present film-substrate 
combination, complete information about the thermal 
constants was not available to calculate R t h from first 
principles, but we are able to show by other means that 
sample heating is not appreciable. Further details of the 
measurement techniques have been published in a previ- 
ous review article and other recent papers pL |(| . All 
measurements were made in zero applied magnetic field 
and the highest self field of the current (~ 300 G) is of 
the order of the lower-critical field H c \ = \%h G M. 



RESULTS AND ANALYSIS 

Fig. |j](a) shows the resistive transitions at different 
electric currents I for the medium sample. The inset 
shows the sample geometry. The horizontal sections of 
the current leads add a small (~ 15 %) series resistance 
to the actual resistance of the bridge. Because j in these 
wide regions is negligible, this resistance freezes out at the 
nominal unshifted T c , making the onset seem to not shift. 
However over the main portions of the curves, there are 
substantial and relatively parallel shifts induced by cur- 
rent. Fig. 0(b) provides a magnified view of the central 
two-thirds portion of the transitions. The dashed line 
represents half the normal-state resistance (i?„) of the 
bridge, which serves as the criterion (resistive-transition- 
midpoint) for defining T c (j). Panels (c) and (d) show 
similar sets of curves for the other two (small and large) 
samples. 

Fig. ^| shows the midpoint T c 's and their correspond- 
ing currents (ranging from 10~ 6 to 10 -2 A) plotted as 
I 2 / 3 (expected for pair-breaking) and as I 2 (expected for 
Joule heating). The shifts are closely proportional to 
I 2 ! 3 rather than to I 2 , showing that heating is not ap- 
preciable (the plots for samples S and L look similar). 
The slope dl 2 ^ 3 /dT c (j) together with Eq. gives a zero- 
temperature depairing current value of 257 mA. Dividing 
this by the cross-sectional area gives a current density of 
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FIG. 1: Resistive transitions of MgBi bridges at different cur- 
rents (values correspond to curves from left to right.). Panels 
(a) and (b) show two windows of the same data. The inset 
in (a) shows the sample geometry and configuration of leads. 
Panels (b), (c), and (d) show the central main portions of 
the transitions for three different sized samples. The right- 
most curves at 1= 1.4 fiA were measured by with a continuous 
DC current; the rest used pulsed signals. The dashed lines 
represents R — R„/2 for each sample. 
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FIG. 2: Shifted transition temperatures at different currents. 
The two y-axes plot the same T c (j) data versus I 2 ^ 3 and I 2 , 
showing adherance to the I 2 ^ 3 law for pair breaking rather 
than the I 2 law for Joule heating. The linear fit (solid line) 
to the I 2/3 plot gives I d (0) = 257 mA (see Eq. 



jd(0) — 2.1 x 10 7 A/cm 2 . The respective values for sam- 
ples S and L are j d (0) = 2.2 x 10 7 and 1.8 x 10 7 A/cm 2 . 
The three values are consistent within the uncertainities 
in the sample dimensions, implying a cross-sectionally 
uniform current density. This is expected for the dissipa- 
tive state of a superconductor (In the fluctuation region 
near the T c (j) boundary and during flux motion-when 
the superconductor is resistive — the current flow becomes 
macroscopically uniform, as in a normal conductor, due 
to the principle of minimum entropy production. This 
has been discussed and verifed elsewhere 0.) and close 
to the T c (j) boundary where A and £ (coherence length) 
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FIG. 3: (a) IV curves for sample M at seven fixed tempera- 
tures (listed for curves going from left to right). Beyond jd, 
the voltage jumps to a linear behavior reflecting the resistance 
of the normal state (slope indicated by the dashed line), (b) 
IV curves for the largest sample in different thermal environ- 
ments to evaluate Joule heating. 

Fig. ||(a) shows current-voltage {IV) characteristics at 
various fixed temperatures for sample M (results for sam- 
ples S and L are similar). As / is increased V remains 
close to zero until some critical value. Above this it shows 
Ohmic behavior V = IR n . Note that at T=35.5K the 
transition is gradual, whereas at the lower temperatures 
it is rather abrupt. This is in part because a type II 
superconducting phase transition changes from second 
order to first order at lower temperatures in the pres- 
ence of a current || and possibly because of a thermal 
component. The "s" shape arises because the external 
circuit feeding the pulsed signal has a source impedance 
R s of about 12 SI. Thus when the sample is driven nor- 
mal, the current will drop discontinuously by the fraction 
R n /(R n + R s ) ~ 20% as observed. Although the amount 
of Joule heating cannot be directly estimated, its sig- 
nificance can be assessed by measuring the IV curves 
in different thermal environments. The previous curves 
in Fig. ||(a) were all measured with the sample in he- 
lium vapour. Such measurements were repeated with 
the samples in superfluid and normal liquid helium, and 
in vacuum, and Fig. ||(b) shows one such set for the 
large sample (its lower surface-to-volume ratio gives it 
the worst thermal resistance). Fig. ||(b) shows no sig- 
nificant systematic influence of the environment on the 
observed value of jd, which would not be the case if Joule 
heating were a serious problem || . 

From such IV characteristics measured at the lowest 
temperature (1.5 K) in superfluid helium, the current 
required to drive the sample normal provides a direct 
lower bound on jd(T « 0) (This lower bound will equal 
jd in the case of uniform current flow. We return to this 
point again later.). For the three samples S, M, and L 
these respective values are jd(0) = 1.9 x 10 7 , 2.0 x 10 7 , 
and 1.7 x 10 7 A/cm 2 , which are consistent with the values 
obtained earlier (2.2 xlO 7 , 2.1 xlO 7 , and 1.8 x 10 7 A/cm 2 ) 
from the shifts in the resistive transitions near T c (Fig. |2j 
and Eq. 

Fig. | shows the values of Id at different temperatures 
obtained from the IV characteristics of Fig. ||(a). Also 
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FIG. 4: Pair-breaking currents from IV curves (Fig. ^) and T c 
shifts (Figs. I] and^j. The solid line represents the theoretical 
curve (Eq. 15) in which ld(0) and T c are fixed by the T c shifts 
(Fig. m; 36 < T < 39,), with no adjustment made over the rest 
of the range (1.5 < T < 36 K). 



order of magnitude to the value of 6.1 x 10 7 A/cm 2 cal- 
culated from Eq. ^ and the published values of -ff c =2500 
G and A = 185 nm from the review on MgB2 by Wang 
et al. j8j, in view of the uncertainities in those param- 
eters. From a technological standpoint, the depairing 
current density of MgB2 is about an order of magnitude 
lower than the high-T c cuprates jy]. The good news is 
that the flux pinning in films is so strong (because of the 
larger coherence length and more isotropically 3-D be- 
havior) that the depinning j c at modest fields appear to 
be within an order of magnitude of jd 12 , whereas for 



the cuprates j c and jd can be separated by two or three 
orders of magnitude |l3j| . 
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shown are the values of Id obtained from the shifts in 
the resistive transition near T c (from Fig. ^). The solid 
line is a plot of Eq. || in which the values of T c (0) and 
id(0) came directly from the observed j 2 / 3 behavior of 
Fig. ||and were not adjusted to fit the other data over the 
extended temperature range, nor was Id(fi) adjusted to 
fit the actual measured value from the IV characteristic 
at low T (i.e., Fig. ||). Nevertheless, the Id(T) data tend 
to follow the general trend of Eq. [| 

An aspect of these data that may seem surprising is 
that even for T <C T c (where the sample width w 3> 
A, £), the average current density in the bridge reaches 
essentially the full jd before the system becomes normal. 
By contrast, previous studies |ic|] of jd in wide low-T c 
type-I bridges found that the current distribution was 
non-uniform and the sample was driven normal when the 
peak j near the surface exceeded jd- Using their model 
with our sample dimensions and parameters, our effective 
jd should have been reduced by a factor of 3, but it is not. 
In our present relatively high-T c type II material, it seems 
that the slight flux motion induced by the self field and 
fluctuations (because of the proximity to the T C (J) phase 
boundary and much higher T) serve to homogenize the 
current distribution. Ironically, such incipient dissipation 
may actually stabilize the flow and permit the average j 
to get closer to jd before the system becomes normal. 

In conclusion, we have measured the fundamental pair- 
breaking critical current density of magnesium diboride 
over the entire temperature range for in-plane current 
transport. The measured jd{T) function is consistent 
with the Ginzburg-Landau form and conforms exactly 
to the AT C cx j 2 / 3 behavior predicted near T c . jd(0) 
obtained from the value of current required to drive 
the sample normal at T — > 0, agrees with the jd{0) 
deduced from the AT C oc j 2 / 3 behavior close to T c . 
The average value for all samples by both methods is 
jd(0) ~ 1.9 ± 0.4 x 10 7 A/cm 2 . This is comparable in 
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